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Abstract 

We show that strongly coupled holographic matter at finite charge density 
can exhibit charge density wave phases which spontaneously break trans- 
lation invariance while preserving time-reversal and parity invariance. We 
show that such phases are possible within Einstein-Maxwell-dilaton the- 
ory in general spacetime dimensions. We also discuss related spatially 
modulated phases when there is an additional coupling to a second vector 
field, possibly with non-zero mass. We discuss how these constructions, 
and others, should be associated with novel spatially modulated ground 
states. 



1 Introduction 



The charge density in a metal is usually highly uniform due to strong Coulombic 
interactions. However, a wide variety of metals exhibit stable charge density wave 
(CDW) phases in which the charge density spontaneously becomes spatially modu- 
lated. These phases were first predicted by Peierls in the context of weakly coupled 
one- dimensional systems [1]. They are also known to occur in strongly correlated 
systems and in some cases the CDWs are present with additional order. For ex- 
ample, in the high temperature cuprate superconductors there are striped phases in 
which CDWs appear with spin density waves, with the latter breaking time-reversal 
invariance (for a review see [2]). In this paper we will discuss CDWs, that preserve 
time-reversal and parity invariance (T and P), within the holographic framework of 
the AdS/CFT correspondence. 

The first constructions of holographic spatially modulated phases, but without 
CDWs, were made in the context of D = 5 Einstein-Maxwell theory with a Chern- 
Simons term [3,4] (see also [5]). At finite charge density and high temperatures, 
the system is described by the standard electrically charged AdS-RN black brane 
corresponding to a spatially homogeneous and isotropic phase. For sufficiently large 
Chern-Simons coupling, the AdS-RN black brane becomes unstable at a critical tem- 
perature and a new branch of black hole solutions appears corresponding, in the dual 
field theory, to a phase with spatially modulated currents with helical order [3]. The 
fully back-reacted black hole solutions were constructed in [4] where it was shown 
that helical current phase is thermodynamically preferred and that the phase transi- 
tion is second order. Furthermore, at zero temperature the black hole solutions of [4] 
revealed new ground states with a helical structure similar to those studied in [6]. 
These helical current phases break P and T. 

In subsequent work, a, D = 4 Einstein-Maxwell theory coupled to a pseudo-scalar 
ip was shown to admit black hole solutions corresponding to phases with spatially 
modulated currents and in addition, CDWs [7]. In this class of models the key 
coupling driving the spatially modulated phase transition is the axion-like coupling 
(fF A F, where F is the field strength of the Maxwell field. The fully back-reacted 
black hole solutions require solving partial differential equations and some results 
have recently appeared in [8]. While the models of [7] realise CDWs the spatially 
modulated currents break P and T. Holographic constructions of spatially modulated 
phases at finite charge density 1 have also been made in [16-23] and P and T are again 

1 Spatially modulated phases preserving the U{\) symmetry have also been shown to exist in the 
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not preserved. 

It is natural to ask, therefore, if breaking of P and/or T is necessary to realise a 
spatially modulated phase at finite charge density in the context of holography. Here 
we will show that it is not 2 . Specifically, we will discuss two general classes of models 
associated with CDW phases, without current density waves, and preserving P and T. 
Our general strategy will be to consider set-ups in which the zero temperature limit 
of the electrically charged black holes that describe the unbroken high temperature 
phase, are domain wall solutions interpolating between some UV fixed point and 
an electrically charged AdS2 X M> D ~ 2 solution in the far IR. We then construct CDW 
type modes that violate the AdSz BF bound; the existence of these modes necessarily 
implies that the finite temperature unbroken-phase black holes become unstable at 
some critical temperature. Furthermore, at this critical temperature the unbroken 
phase black holes will admit a static normalisable mode corresponding to the existence 
of a new branch of electrically charged black holes dual to the CDW phase. 

The first class of models, discussed in section 2, involve Einstein-Maxwell theory 
coupled to a single scalar field (f>, often called the dilaton. These models have a single 
conserved U(l) symmetry in the dual field theory. We will see that the relevant CDW 
mode involves the gauge-field, the dilaton and the metric. 

The second class of models, discussed in section 3, involves Einstein-Maxwell- 
dilaton models with an additional coupling to a second vector field which can either 
be massive or massless. This class of models is somewhat simpler than the first in 
that the standard electrically charged AdS-RN black brane solution solves the equa- 
tions of motion. The AdS-RN black branes describe the high temperature unbroken 
phase of a dual CFT when held at finite charge density with respect to the U(l) 
symmetry associated with the Maxwell field. When the second vector field is mass- 
less the dual theory has a second global U(l) symmetry. We show that there can 
be spatially modulated modes just involving the dilaton and the second vector field, 
which correspond, in the massless case, to CDWs for the second U(l) symmetry. 

In section 4 we conclude with some discussion on how our constructions should be 
associated with novel spatially modulated ground states. One route is to follow the 
spatially modulated phases down to zero temperature. A different route, following 
[25], utilises AdS2 solutions which contain spatially modulated modes that are dual 
to relevant operators in the RG sense. 

presence of magnetic fields in [9—15] , where the magnetic field itself breaks T. 
2 For a discussion of constructions exploiting magnetic gaugings, see [24]. 
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2 Einstein-Maxwell-Dilaton models 



In this section we consider Einstein-Maxwell-dilaton models in D spacetime dimen- 
sions, which couple gravity to a gauge-field A, with field strength F = dA, and a 
scalar field 0, the "dilaton" . The Lagrangian density is given by 

£ = ^_ y (0)-I(d0) 2 -|r(0)F 2 , (2.1) 

where F 2 = F^F^ and V, t are functions that we will partially restrict below. The 
equations of motion derived from the action (2.1) can be written as 

Rfxu 9/JiV ^^^(f) d u (f) (^t L P^ u ^ 2(D— 2) 9 ^v^pa-F^ j , 

V M (rFH = 0, 

V 2 - V - \r'F 2 = . (2.2) 

This class of models has been studied in a holographic context in [26] (see also e.g. 

[27-35]). 

We will consider particular classes of models that admit electrically charged 
AdSi x ]R D ~ 2 solutions of the form 

ds 2 = L 2 (ds 2 (AdS 2 ) + dx\ H h dx 2 D _ 2 ) , 

F = EVo\(AdS 2 ) , 

= 0o, (2.3) 

where E, O , L are constants and we have written the metric and field strength using 
two dimensional anti de-Sitter space with unit radius. Using the equations of motion 
(2.2) this implies the conditions 

V{cp ) t> ) = -t(¥>o) V'icpo), 

V(<p )<0, r(^ )>0, (2.4) 

and 

V{<po) -V 0Po)t(<Po) 

We will take E > 0. We will need the expansion of the potential V and the function 
r around the value (fo. In general, it can be brought to the convenient form 

V = V {1 - Ti (<p - Ifo) - — {if - (fo) +■■■), 
T = T (l+T 1 (if- (fo) - y (<p - tp ) 2 H J , 

v <0, r >0, (2.6) 
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where we have incorporated the conditions in (2.4). 

We now wish to study linearized perturbations around the background (2.3). For 
simplicity we will continue with D = 4 but we will also quote the final key result for 
D = 5. More specifically, using a coordinate system for AdS 2 such that 

dr 2 

we are interested in the perturbation 



ds 2 (AdS 2 ) = -r 2 dt 2 + — , (2.7) 



5g tt = L 2 r 2 e~^h tt (r) cos(k Xl ) , 
5g XiXi = L 2 e~ lu)t h x . x . (r) cos{kx{) , 
6g tXl = L 2 e~ lU)t h tXl (r) sva.{kx x ) , 
5 A t = -Ee- iut a t (r) cos(k Xl ) , 
SA Xl = -Ee~ iut a Xl (r) wa.{kx x ) , 

= e - iuJt h(r) cos(k Xl ) , (2.8) 

with i — 1,2, which involves seven functions of the radius {h t t,h XiXi ,ht Xl ,a t ,a Xl ,h} 
in a self consistent manner, as we shall see. 

The linearised perturbation is in a radial gauge but there is still some remaining 
gauge freedom. Specifically, the combined co-ordinate transformation — > + 5x^ 
and U(l) gauge transformation A^ — > A^ + d^5A with 

dr — r f , 
5xi = - \nrd Xl f + g 2 , 
Sx 2 = , 

5A = --d t f + Eg 3 , (2.9) 
r 

keeps the form of the perturbation (2.8) invariant when / = f(t, xi) and gi = gi(t, x\). 
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Indeed for such gauge transformations we have 

6g XlXl = 2L 2 (~\iard 2 x J + d Xl g 2 ) , 

$9x2X2 = ) 

8g txi = L 2 (^d t d Xl f - r 2 d xi g x - \nrd t d Xl f + d t g 2 

SA t = E f-rf - \ l d 2 f - rd t9l + d t g 3 
\ 2r 

5A X1 =E (-±-d t d x J + d Xl g 3 -rd Xigi ) ; (2.10) 

and we see that two gauge invariant combinations are given by 

1 2 

jY^Sgu - — d r 5A t , 5g X2X2 . (2.11) 

This will be useful shortly. 

We now return to the linearised perturbation (2.8) about the AdSz x solution 
(2.3). After substituting into (2.2) we obtain three ordinary differential equations 
from the gauge field equation of motion, seven from Einstein's equations and one 
from the scalar equation of motion: 

- 2r 2 a' t + r 3 (g' XiXi + g' X2X2 ) + u 2 {g XlXl + g X2X2 ) - 2ikug tXl - r 4 g" t 

- 3r 3 flk + k 2 r 2 g tt - r 2 g tt = , 

- iuj (-r (g' xixi + g' X2X2 ) + g xixi + g X2X2 ) + krg' txi - 2kg txi = , 
2a' t - r (r (g£ ixi + g" x%X2 - gQ + g' XlXl + g' X2X2 - 3g' tt ) + g tt = , 

r 2 (2a' t + 2rg' XiXi + 2r x h + r 2 g" iXi + (1 + k 2 )g tt - k 2 g X2X2 ) + u 2 g XlXl 

- 2ikug tXl = , 
9x2x2 + r 2 (2a' t + 2rg' X2X2 +2r x h- k 2 g X2X2 + r 2 g X2X2 + g a ) = , 
r 2 {2T X a' t + r 2 h" + 2rti - k 2 h - (r 2 + v 2 )h + T X g tt ) + co 2 h = , 
iuj (2a' t + g XlXl + g X2X2 + g tt + 2r x h) + 2kr 2 a' Xi + 2kg txi = , 
2ikua Xl - 2k 2 a t + r 2 (2< + g' xixi + g' X2X2 + g' tt + 2 n ti) = , 
r 2 { 2a ' Xl + O'U) + ikug X2X2 = , 

- k (2a t + r(r (g' tt - g' X2X2 ) + g u )) + iu (2a Xl + g' tXi ) = , 
(r 2 a Xl + 2ra' Xl + g' tXl ) + u 2 a Xl + ikua t = . (2-12) 



r 
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Some analysis now shows that four of these ODEs are implied by the remaining seven, 
as expected from the Bianchi identity for the Einstein tensor. Furthermore, we find 
that in this system of seven ODEs the three functions {h, a t , h X2X2 } enter with second 
order derivatives in r while the four functions {h xixi ,h txi ,h tt ,a xi } enter with first 
derivatives. This indicates that there are three propagating degrees of freedom and 
four constraints. 

To proceed we now exploit our previous discussion on gauge invariance and intro- 
duce the variables 

<5> l = h tt + 2a' t , <5>2 = g X2X2 , &3 = h, (2.13) 

and find that the system of seven ODEs can be written in terms of {$1, $ 2 , $3} and 
{K 1X1 , h txi , h tt , a xi } in the form: 

w 2 $ 3 + r 3 (r$ 3 ' + 2$' 3 ) + r 2 (n$i - (k 2 + v 2 + r 2 )$ 3 ) = , 

w 2 $ 2 + r 3 (r$' 2 ' + 2$' 2 ) + r 2 - k 2 <& 2 + 2<n$ 3 ) = , 

2iku(ra Xl + g tXl ) + 2k 2 r 2 a' t - 2k 2 ra t - co 2 g XlXl + 2r 3 n$ 3 - r 2 (2n$ 3 + $1) 

- fcV$i + k 2 r 2 § 2 + r 3 ^ - w 2 $ 2 = , 

- 2k 2 r 2 a' t + r 3 g' xixi + u 2 g xixi - 2ikuug txi + 2r 2 r 1 $ 3 + (k 2 + l) r 2 ^ 

+ (cu 2 - k 2 r 2 ) $ 2 + r 3 $ 2 = , 
iu (-2k 2 r 2 a' t + 2r 2 ri$ 3 + (k 2 + l) r 2 $ x - fc 2 r 2 $ 2 + r 2 $ 2 + w 2 $ 2 + (r 2 + w 2 ) g XlXl ) 

-kr 3 g' txi +2k (r 2 + u 2 )g tXl =0, 
2ik (k 2 + 1) r 3 coa Xl + 2k 2 r 5 a" + 2A; 2 r 4 (A; 2 + l)a' t + 2k 2 rWa' t 

- 2k 2 (k 2 + 1) r 3 a t - u 2 ((k 2 + 1) r 2 + w 2 ) g xixi + 2^0^ ((/c 2 + l) r 2 + u 2 ) 
+ 2A; 2 r 5 r 1 $ 3 - 2r 2 r x (k 2 r 2 + u 2 ) $ 3 + (k 4 r 4 - u 2 (r 2 + u 2 )) $ 2 + k 2 r 5 & 2 

-r 2 (k 2 ((k 2 + 2) t 2 + uj 2 )+uj 2 ) ^ = , 
iw^x!*! + 2n$ 3 + $! + $ 2 ) + 2A; (r 2 a^ + = . (2.14) 

The first two equations comprise part of a second order system for the three gauge- 
invariant fields {$1, $ 2 , $ 3 }. We can solve the third equation for Differentiating 
this expression and then using the remaining equations, we find the remaining second 
order equation. Introducing the three-vector V T = ($1, $ 2 ,$ 3 ) we can write these 
in matrix form as the standard system of three mixed modes propagating on AdS 2 '- 

(j 2 - + r 2 d 2 .+2rd r ^jY-M 2 Y = 0, (2.15) 
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with mass matrix 



M 



( 2 + r 2 + k 2 -2k 2 2-n (2 - k 2 - r 2 - v 2 ) \ 

-1 k 2 -2r x 

\ -ri k 2 + v 2 + r 2 J 



(2.16) 



In addition, once we have specified these three fields the remaining four functions 
are determined by solving first order equations arising from (2.14) and demanding 
regularity at the Poincare horizon of the AdS 2 space at r = 0. 

To summarise, for D = 4 we have shown that the linearised perturbation (2.8) 
corresponds to three propagating modes in AdS 2 with mass matrix (2.16). We note 
that the mass matrix (2.16) only depends on two parameters appearing in (2.6), T\ 
and r 2 + v 2 . As a check we note that when t% = the scalar field decouples from the 
metric and gauge-field fluctuations and we obtain the mass-squared eigenvalues 

ml = k 2 + v 2 + t 2 , 



m 



± 



1 + k 2 ±Vl + 2k 2 . (2.17) 



In particular, the masses m\ agree with previous studies of longitudinal perturbations 
of Einstein-Maxwell theory around its AdS 2 x IR 2 solution [36]. 

The AdS 2 BF bound is violated if any of the three eigenvalues m 2 of the mass 
matrix (2.16) is such that m 2 < —1/4. The signal for spatially modulated phases are 
modes that violate the BF bound with the smallest mass-squared occurring at k ^ 0. 
This is easily achieved by suitable choices of the parameters T\ and r 2 + v 2 , including 
cases where there BF violating modes only occur for k ^ 0. Indeed, the mass squared 
eigenvalues when k = are simply 

ml = 0, m 2 = 2, m 2 3 = 2t 2 + r 2 + v 2 . (2.18) 

Demanding that 2t 2 + t 2 + v 2 > —1/4 there are still broad choices for the parameters 
which have modes with fc ^ that violate the BF bound. Note from (2.17) that 
while there can be instabilities when n = they will not be spatially modulated. 
To illustrate we can consider the specific choice of V, t given by 

V = v e-~" t> , r = e^, (2.19) 

where 7 is a constant, that were studied in [26]. The eigenvalues of the mass matrix 
are given by k 2 , 1 + k 2 ± a/1 + 2(1 + 7 2 )A; 2 . At k = eigenvalues are 0, and 2 and 
hence do not violate the BF bound. If 7 > 1 there are BF violating modes for a range 
of k 7^ 0, with the minimum mass-squared occurring at k 2 = 7 2 (2 + 7 2 )/(2(l + 7 2 )). 
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We emphasise that the spatially modulated BF violating modes that we have 
identified are associated with CDWs preserving P and T. To see this 3 we first note 
that the mode that is relevant for this discussion is static with u> = 0. From (2.14) we 
deduce that these modes have h txi = a xi = and hence the linearised perturbation 
in (2.8) has Sg txi = SA Xl = 0. Thus, the gauge-field part of the linearised perturba- 
tion only involves a spatially modulated A t component corresponding to a spatially 
modulated CDW in the dual theory. It is also clear from the perturbation in (2.8) 
that the CDW preserves P and T. Indeed it is worth noting that the mass-matrix 
in (2.8) only depends on k 2 in contrast to other cases that have been analysed in the 
literature (e.g. equation (3.17) of [3] and equation (2.17) of [ ]) where it depends 
linearly on k. 

Finally, we have focussed on Einstein-Maxwell-dilaton theory in D = 4 spacetime 
dimensions, but very similar comments apply to other D. For example, when D = 5, 
the analogue of the mass matrix given in (2.16) is presented in appendix A. 

3 Einstein-Maxwell-dilaton-vector models 

We now consider a theory in D bulk space-time dimensions which couples gravity to 
a massless U(l) gauge field A, a scalar dilaton field <fi and an additional vector field 
B of mass m v . The Lagrangian density is given by 

C=R- \(d<P) 2 - V (0) - \t (0) F 2 - \v(<P)G 2 - \m 2 v B 2 - \u (0) FG , (3.1) 

where F = dA and G = dB and FG = F^G^ . The corresponding equations of 
motion are given by 

R p» = lh V 9^u + \m 2 v B^B u + \d^dA + \t (F^F/ - ^L^g^F^F^ 

-\- yG upG^ 2(D—2) 9 [ivG p a G^ j -\- u i^G ^jj 1 F v ^p 2(0—2) 9 ^uG p a F^ ^ , 
Vp {tF^ + uG^) = 0, 

(vG^ + u F» v ) - m 2 v B v = , 
V 2 -V- \t'F 2 - \v'G 2 - \u'FG = . (3.2) 

3 We will elaborate on this discussion in a simpler setting in the next section in the context of 
another class of models. 
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We will assume that the functions V, t, u and v have the following expansion 

t(0) = l-^L 2 2 + --- , 
u(<t>) = ~^ sL ^ + ' ' ' > 

v{<j)) = l + --- . (3.3) 

Then AdSp, with radius squared £ 2 , where 

£ 2 = L 2 (D- 1) (D - 2) , (3.4) 

and = A = B = 0, solves the equations of motion and is dual to a CFT in 
D — 1 space-time dimensions. The gauge-field A is dual to a conserved current for a 
global [/ (1) symmetry, while <ft and B are dual to neutral scalar and vector operators, 
respectively. In the special case that m v = the dual CFT has a second global U(l) 
symmetry and B is dual to the conserved current. 

We are interested in the CFT held at fixed chemical potential fi with respect to 
the global U(l) symmetry associated with the gauge-field A. The high temperature 
phase is described by the standard electrically charged AdS-RN black brane solution 
with = B = 0. At zero temperature, the solution interpolates between AdSo in 
the UV and AdS2 x M, D ~ 2 in the IR, with the latter solution given by 

ds\ = L 2 (ds 2 (AdS 2 ) + dx\ H h dx 2 D _ 2 ) , 

F = Evo\ (AdS 2 ) , 

<j) = G = 0. (3.5) 

where E = \/2L. 

In the following subsection we will first examine CDW type instabilities of this 
AdS 2 x M. D ~ 2 solution. Such instabilities imply that electrically charged AdS-RN will 
have analogous instabilities at finite temperature, and these will be analysed in the 
subsequent subsection. We also note that the spatially modulated instabilities that 
we saw in the last section are not present in the models that we consider in this 
section because in (3.3) we have assumed t'(0) = 0. 
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3.1 Instabilities of the AdS2 X R D 2 solution 

We consider the following linearised perturbation about (3.5), using the coordinates 
for AdS 2 given in (2.7), 

56 = e -^+ ik ^ $( r ) , 



<)D = f ( r '%(r)dt + iub r {r)% ) . (3.6) 



We now substitute into the equations of motion (3.2). The x\ component of the 
equation of motion for the gauge-field B implies the constraint 4 

kcu(b t + b' r ) = . (3.7) 

We also obtain a coupled system of equations which we can write in matrix form 



to 2 

+ r 2 d 2 + 2rd r ) V -L 2 M 2 Y = 0, (3.S 



where V T = ($, b r ) and the mass matrix is given by 



M 2 = i ■n 2 + n + s 2 +p 2 -s{m 2 v +p 2 ) 
— s ml+p 2 



with p = k/L. The matrix (3.9) yields the AdS2 mass spectrum 



2 1 / ~ 2 i 2 i 2\ , 2 
m ± = 2 (. m s + m « + S ) + P 



±^\/(^-^) 2 + 2 (^ + ^ + 2p 2 ) S 2 + S 4 , (3.10) 

where 

m 2 = m 2 + n . (3.11) 

It is now straightforward to choose parameters such that there are modes violating 
the BF bound associated with spatially modulated phases. As a specific example, 
which we will return to in the next subsection, if we set D = 4, L 2 = 1/24, m 2 = —8, 
ml = 0, n = —96, s = 16.2 we find L 2 m 2 _ < —1/4 in the range 

3-52 < |p| < 8.92. (3.12) 



4 Note that when ^ the equations of motion (3.2) imply that W^SB^ = and hence 
co(b t + b' r ) = 0. We also note that when ml = we can work in a gauge with 8B r = if desired. 
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Note that when u = the spatially modulated BF violating modes involves the 
scalar field and only the time component of the gauge-field, B t . To be more explicit, 
consider the static mode 

6(f) = e ikxi v ir x , 

5B = e lkxi v 2 r x+1 dt, (3.13) 

where Vi,v% are constants and A corresponds to the scaling dimension of an operator 
in the one- dimensional CFT dual to the AdS 2 x R D ~ 2 solution. We obtain a solution 
to the linearised equations of motion provided that 




This can be solved for A and it is straightforward to check that A(A + 1) = m| as 
expected. As usual, the BF bound is violated when A becomes imaginary, the onset 
of which happens when A = —1/2. These BF violating modes are associated with 
spatially modulated phases that preserve P and T. Furthermore, in the special case 
that m 2 = 0, when there is a second global U(l) symmetry in the dual CFT, the 
unstable modes correspond to CDWs for the second Z7 (1). We will make this more 
explicit in the next subsection. 

3.2 AdS-RN black hole instabilities 

In this subsection we will explicitly show, for an illustrative case, that the instabilities 
that we deduced for the AdS 2 x M D ~ 2 solutions are associated with instabilities of the 
finite temperature AdS-RN black brane solution. In particular, we will calculate the 
value of the critical temperature at which the AdS-RN black brane becomes unstable 
as a function of wave number, finding the usual "bell-curve" type behaviour. 

We will work in D = 4. The electrically charged AdS-RN black brane is given by 

ds 2 = -fdt 2 + ^y- + r 2 {dx\ + dx 2 2 ) , 

A = (1 - T -±)dt, 
r 

<j) = B = 0, (3.15) 

with 

f = r *( r A + ±)r-± + ±. (3.16) 



£ 2 \£ 2 4 J r 4r 2 
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This describes the high temperature phase of the dual d = 3 CFT at finite tempera- 
ture T and non-vanishing chemical potential fi with respect to the global symmetry 
corresponding to the gauge-field A. Notice that we have chosen to scale so that \i = 1 
and that T = (12r 2 - e 2 )/(16irr + £ 2 ). We will also set m 2 v = so that the CFT has 
a second global symmetry, corresponding to the gauge-field B. The instabilities that 
we discuss correspond to CDW phases associated with this second U(l) symmetry. 
Specifically, we consider the following linearised perturbation 

5(f) = 4>{r) cos(fcxi) , 
SB = b t {r) cos(kx 1 )dt . (3.17) 

We will take the scalar field to have mass given by m 2 = —2/£ 2 . This corresponds 
to having an operator in the dual d = 3 CFT with scaling dimension A = 1, 2 and 
we will choose the boundary conditions so that A = 2. Since we are interested in 
instabilities associated with phases that spontaneously break translation invariance, 
we consider the following UV expansion as r — > oo, in which the sources are set to 
zero: 

^(r)«^ + ---, &t(r) «- + '••■ (3.18) 

We will also demand regularity on the black hole horizon at r = r + , by demanding 
that the functions admit the analytic expansion: 

<!>(r) « 0(o) + O (r - r+) , 

h{r) « 6(o) (r -r+)+0{r- r + f . (3.19) 

The linearised equations of motion lead to two second order differential equations 
for (f>(r), b t (r) and so a solution is specified by four integration constants. For a given 
k we have five parameters, r + , fa, q, <f>(p), 6( ) entering the ODEs. Since the ODEs 
are linear, we can always scale one of the parameters to unity. This means that 
for a given k, we expect solutions to exist, if at all, for specific critical values of T. 
These are precisely the static modes at the onset of the instability. In figure 1 we 
have plotted the critical temperature versus wave number for these modes for the 
specific choice of parameters given by n = —96, s = 16.2 and I = 1/2. Recall that 
we used the same parameters in the last subsection and we also note that in contrast 
to the last subsection, the equations now depend on m 2 and n individually and not 
just in the combination (3.11). To make a comparison with the wave-numbers of the 
AdS2 x 1R 2 analysis in the last subsection given in (3.12), we should rescale the spatial 
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Figure 1: A plot of the critical temperature versus wave number at which the 
electrically-charged AdS-RN black brane becomes unstable to the formation of a 
CDW phase preserving P and T. The plot is for the model (3.1), (3.3) with n = —96, 
s = 16.2, m 2 s = —8 and ml = 0. For these values the highest critical temperature 
occurs at T c //i « 0.011 with \k c \//i « 0.645. 

coordinates by a factor of \/2 leading to dividing p by 4\/3 so that (3.12) corresponds 
to 0.508 < | A; | < 1.29, in good agreement with figure 1. 

The static mode with the highest critical temperature, T c , associated with wave- 
number k c , corresponds to the onset of CDW phases in the dual CFT. Indeed, for 

( B) 

temperatures just below T c the charge density j\ for the second Z7 (1), corresponding 
to the gauge-field B, becomes spatially modulated with 



At the same time the operator O^, dual to the scalar field 0, also spontaneously 
acquires an expectation value given by ~ (f>2 cos(/c c Xi). It is clear that this new 
CDW phase preserves P and T. 

4 Discussion 

We have shown that two classes of holographic models with electrically charged 
AdS2 x ]R D_2 solutions can have spatially modulated BF violating modes that are 
associated with CDWs. Specifically, these modes demonstrate the existence of rich 
classes of spatially modulated black hole solutions that are dual to the appearance of 
CDW phases, preserving both P and T, in holographic matter at finite charge density. 
For the models discussed in section 3 with two £7(1) vector fields, we deduced the 



(jt ) ~ gcos(fc c xi) . 



(3.20) 
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critical temperature at which these black hole solutions will appear. Similar calcula- 
tions are also possible for the Einstein-Maxwell-dilaton theory considered in section 
2, by choosing specific functions V, t, constructing the unbroken phase black holes 
and then analysing the linearised perturbations. 

Constructing the fully back reacted black hole solutions will require solving non- 
linear partial differential equations. The simplest black hole solutions will depend 
on two variables, the radius r and the coordinate xi, being static and translationally 
invariant in the remaining spatial coordinates. However, since at the linearised level 
we can superimpose the static normalisable modes, there will also be more elaborate 
black hole solutions that depend on all of the spatial coordinates. Only a detailed 
analysis will reveal which periodic structure is thermodynamically preferred. It would 
be particularly interesting to construct the black hole solutions all the way down to 
zero temperature so that ground states of the system can be identified. While it is 
possible that the spatial modulation disappears at zero temperature it seems much 
more likely that the generic ground states will be spatially modulated CDWs. 

The results of this paper, combined with those of [25], also suggest other con- 
structions of novel spatially modulated holographic ground states. We suppose the 
AdS 2 x R D ~ 2 solution still arises as the IR limit of a domain wall solution inter- 
polating from some holographic behaviour in the UV. But now we assume that the 
AdS 2 x IR- 0-2 solution is stable, without any modes violating the BF bound. In ad- 
dition we assume 5 that the AdS 2 x M. D ~ 2 IR fixed point solution does not have any 
relevant operators with k — 0, but does have relevant operators with k ^ 0. If we 
now consider deforming the UV fixed point of the domain wall by a spatially homo- 
geneous and isotropic deformation, the IR fixed point will be stable under RG flow, 
because of the absence of k = relevant modes in the IR. However, if we switch on a 
suitable spatially modulated deformation in the UV (such as a spatially modulated 
chemical potential) the RG flow will be destabilised due to the relevant IR operators 
with k 7^ 0. It seems likely that such constructions will also lead to rich classes of 
spatially modulated ground states. As in [25], these constructions should include 
metal-insulator transitions where the metallic phase is described by the AdS 2 x R D ~ 2 
IR behaviour and the insulating phase by the putative spatially modulated ground 
state. The spatial modulation of the charge density suggests that such insulating 
phases can be interpreted as holographic Mott insulators. 

The holographic CDW phases that we have found in this paper were all in the 
context of models in which the unbroken phase black holes have a zero temperature 

5 A concrete example is provide by the model (2.19) with < 7 < 1. 
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limit containing an AdS 2 factor in the IR. This was purely a technical simplification 
and similar spatially modulated phases should also occur in many other situations. 
For example, within Einstein-Maxwell-dilaton theory there can be black holes whose 
zero temperature limit interpolates between some UV behaviour and a hyperscaling 
violating ground state [26, 33, 34] in the far IR and these should manifest similar 
phases, leading to even more spatially modulated ground states at zero temperature. 
Support for the conjecture made in [21] that the generic holographic states at finite 
density and/or in a magnetic field are spatially modulated continues to accumulate. 
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A D = 5 mass matrix 

For Einstein-Maxwell-dilaton theory in D space-time dimensions the conditions (2.4)- 
(2.6) are still applicable. The relevant perturbation about the AdS 2 x IR 3 solution 
is 



5g tt = L 2 r 2 e^htt (r) cos(fczi) 
5g x . x . = L 2 e~ lU)t h x . x . (r) cos(fcxi) 
8g txi = L 2 e~ iujt h txi (r) sin(fcxi) , 

5 A t = —Ee~ iuJt a t (r) cos(A;a;i) , 
SA Xl = —Ee~ iuJt a Xl (r) sin(fcxi) . 



-iwt 



h (r) cos (k xi) 



(A.1) 



with h X2X2 = h X3X:i . The analysis is very similar to that of the D = 4 case. Using 
the same definition of the three scalar AdS% fields given in (2.13) we obtain the mass 
matrix 



M 



( 2 + 2 t 2 + k 2 -Ak 2 2n (2 - k 2 - r 2 - v 2 ) \ 
-2/3 k 2 -4rx/3 

\ -n o k 2 + v 2 + r 2 J 



(A.2) 
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